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The congruent transformation of the electronic Hamiltonian is developed to address the electron 
correlation problem in many-electron systems. The central strategy presented in this method is 
to perform transformation on the electronic Hamiltonian for approximate removal of the Coulomb 
singularity. The principle difference between the present method and the transcorrelated method 
of Handy and Boys is that the congruent transformation preserves the Hermitian property of the 
Hamiltonian. The congruent transformation is carried out using explicitly correlated functions and 
the optimum correlated transform Hamiltonian is obtained by performing a search over a set of 
transformation functions. The ansatz of the transformation functions are selected to facilitate ana- 
lytical evaluation of all the resulting integrals. The ground state energy is obtained variationally by 
performing a full configuration interaction (FCI) calculation on the congruent transformed Hamil- 
tonian. Computed results on well-studied benchmark systems show that for the identical basis 
functions, the energy from the congruent transformed Hamiltonian is significantly lower than the 
conventional FCI calculations. Since the number of configuration state functions in the FCI calcula- 
tion increases rapidly with the size of the 1-particle basis set, the results indicate that the congruent 
transformed Hamiltonian provides a viable alternative to obtain FCI quality energy using a smaller 
underlying 1-particle basis set. 
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I. INTRODUCTION 

The form of the many-electron wavefunction in the 
proximity of the electron-electron and electron-nuclear 
coalescence point plays a critical role in accurate deter- 
mination of the ground and excited state energies. Al- 
though, the precise structure of the many-electron wave- 
function continues to be elusive, the form of the exact 
wavefunction at the coalescence point is well understood 
and is given by the Kato cusp condition. [1-4] In the 
many-electron wavefunction, the electron-nuclear cusp 
condition can be incorporated by using Slater-type or- 
bitals (STOs). For calculations involving Gaussian- type 
orbitals (GTOs), the one-electron basis can be improved 
iteratively by adding GTOs with increasing angular mo- 
mentum quantum number. [ ] The subject of convergence 
of single-particle basis has been analyzed extensively us- 
ing both analytical and numerical techniques. [6—8] The 
electron-electron cusp has been the focus of intense re- 
search because of its direct relation to the electron corre- 
lation problem and accurate description of the Coulomb 
and Fermi hole. [9-13] However, unlike the electron- 
nuclear cusp, atom-centered basis functions are not ideal 
for accurate description of the many-electron wavefunc- 
tion near the electron-electron cusp. [14-16] Indeed it 
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has been shown that the slow convergence of a full con- 
figuration interaction (FCI) calculation with respect to 
the one-particle basis is related to the inadequate treat- 
ment of the electron-electron cusp. [15] The solution is 
to include explicit r 12 dependence in the form of the 
wavefunction, and there is a large assortment of quan- 
tum chemical methods that have incorporated this ap- 
proach. For example, in the variational Monte Carlo 
(VMC) method, the Jastrow function is used for includ- 
ing explicit r\2 terms in the trial wavefunction. [2, 4] 
The form of the Jastrow is chosen to ensure that the 
electron-electron and electron-nuclear Kato cusp condi- 
tions are satisfied. The parameters in the Jastrow func- 
tion are obtained by minimizing the linear combination 
of energy and its variance. Because of the complicated 
mathematical form of the Jastrow function it is not pos- 
sible to evaluate the integral over the electronic coordi- 
nates analytically, and a stochastic numerical method is 
used for computation of the energy. Recently, Morales et 
al. performed highly accurate multi-determinant VMC 
calculations on water. [ I 7] A detailed review of various 
applications of Quantum Monte Carlo (QMC) methods 
in physics and chemistry can be found in Ref. 18-20. 
Explicitly correlated methods have also been developed 
for post Hartree-Fock schemes such as perturbation the- 
ory (MP2-R12), coupled-cluster methods (CC-R12), and 
multirefcrence CI schemes (R12-MRCI). These methods 
introduce the electron-hole interparticle distance directly 
into the calculation in order to increase the accuracy of 
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the calculations. The field of explicitly correlated meth- 
ods for electronic structure calculation has been reviewed 
and a detailed description of various methods can be 
found in Ref. 14, 15, and 21. A common feature of 
the R12 and F12 methods discussed above is that they 
all involve analytical computation of the T\i correlation 
function. Recently, Chinnamsetty and coworkers have 
presented an interesting study that compared and con- 
trasted QMC with various F12 methods. [22] 

A different strategy known as the transcorrelated 
method was developed by Handy and Boys in 1969. [23] 
The basic idea of the transcorrelated method is to remove 
the electron-electron Coulomb singularity by performing 
similarity transformation on the Hamiltonian using an 
explicitly correlated function. The method was later ex- 
tended by Ten-no to treat the electron-electron cusp us- 
ing Guassian geminal functions and was applied to chem- 
ical systems. [8, 24] The transcorrelated method has also 
been combined with other methods such as QMC [25] 
and coupled-cluster theory [26] and has been used to 
study electron correlation in periodic systems. [ 7 ] One 
of the defining characteristics of this method is that the 
transcorrelated Hamiltonian is not Hermitian and there- 
fore is not required to be bounded from below by the 
exact ground state energy. The correlation function can 
be obtained either by minimizing the energy variance of 
the transcorrelated Hamiltonian [28, 29] or by requiring 
the correlation function to satisfy the electron-electron 
cusp condition. 

The focus of the present work is to address the non- 
Hermitian property of the transcorrelated Hamiltonian 
by replacing the similarity transformation by congruent 
transformation. [30-32] By performing congruent trans- 
formation, we preserve the Hermitian property of the 
electronic Hamiltonian which allows us to use a stan- 
dard electronic structure method such as configuration 
interaction method to minimize the total energy. The 
remainder of the paper describes the theoretical develop- 
ment and the implementation details of the method. The 
derivation of the congruent transformed Hamiltonian is 
presented in Sec. II. Details of performing FCI calcula- 
tions using the congruent transformed Hamiltonian and 
interfacing it with existing FCI methods are presented in 
Sec. II A and II B. Benchmark calculations using the con- 
gruent transformed Hamiltonian are presented in Sec. III. 
The analysis of the results and conclusions are presented 
in Sec. IV. 



II. CONGRUENT TRANSFORMED 
HAMILTONIAN 

The congruent transformed (CT) Hamiltonian H is de- 
fined by performing the following transformation [30-32] 



nian with respect to any trial wavefunction is given as 



H = G^HG, 



(1) 



where G is an explicitly correlated function which will be 
defined later. The expectation value of the CT Hamilto- 
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where 1 = G'IG. The above expression is mathemat- 
ically equivalent to calculating the expectation value of 
the electronic Hamiltonian using a correlated wavefunc- 
tion and is bounded from below by the exact ground state 
energy -E C xact < E^. The optimized energy associated 
with the CT Hamiltonian is obtained by performing a 
minimization with respect to the trial wavefunction and 
explicitly correlated function, 



CT 



min min E^I^t, G]. 

* T G 



(3) 



The optimization of the correlation function G and the 
trial wavefunction is conducted in two steps. In the 
first step, the form of the trial function is kept fixed to a 
single Slater determinant and the parameters of the gem- 
inal functions are determined by minimizing the geminal 
parameters and the molecular orbitals. In the second 
step, the minimized geminal function G m i n is kept fixed 
and the trial wavefunction V&t is minimized. The steps 
involved are described by the following equation 



E[G min ] = min £ t [$ S d,G], 

G,* S d 

E CT = min£ , T [*T,G min ]. 



(4) 
(5) 



The optimization of the correlation function and the trial 
wavefunction are described in the following subsections. 



A. Optimization of the correlation function 

The choice of the correlation function G plays an im- 
portant part in the implementation of the method for 
practical applications. In principle, a variety of corre- 
lated functions such as two and three-body Jastrow func- 
tions can be used. However, the matrix elements associ- 
ated with these functions cannot be integrated analyti- 
cally and one has to use numerical techniques such as the 
VMC method to calculate the integrals. In the present 
work, Gaussian-type geminal (GTG) functions were used 
for the correlated functions. The GTG functions were in- 
troduced by Boys [33, 34] and Singer [35], and have been 
used extensively in explicitly correlated methods. [36- 
41] Slater determinants augmented with GTG functions 
have been used to study electron-electron and electron- 
proton systems. The integrals involving GTG functions 
with GTOs can be performed analytically and have been 
derived earlier. [33, 34, 42, 43] The form of the correlated 
function used in the following calculations is defined as 
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k = l 



where N is number of electrons and N„ is the number of 
Gaussian functions. The geminal coefficients {b k ,j k } in 
the GTG function are determined variationally. In the 
limit of G — > 1, the energy i?[G m i n ] becomes equal to the 
Hartree-Fock energy. 



£?hf = lim E[G min ], 

G->1 



(8) 



As a consequence, the HF energy is the upper bound to 
the geminal minimization process 



E[G m i n ] < Ehf- 



(9) 



The transformed Hamiltonian is expanded as a sum of 
2-6 particle operators as shown below 

^ = EEEfK™)M%(M) (io) 

i<j k m<n 
i<j k<l m<n 

= 2 + 3 + 0± + 5 + 6 (11) 

where, the operators {O n , n — 2, . . . , 6} are defined by 
collecting all two, three, four, five and six particle opera- 
tors obtained by expanding the summation in Eq. (10). 
Specifically, 



2 = Y,h 2 (i,j), 

03 = ^2 h 3 (i,j,k), 

i<j<k 

4 = ^2 hi{i,j,k,l), 

i<j<k<l 

5 = ^2 h 5(i,j,k,l,m), 



(12) 
(13) 
(14) 
(15) 



i<j<fc<2<m 

6 = ^2 h 6 (i,j,k,l,m,n). (16) 

i<j <k<l<m<n 

The exact form of the operators {h n , n = 2, . . . , 6} have 
been derived earlier and are not duplicated here. [44] 
It should be emphasized that the operators {h n ,n — 
2, ... ,6} are defined so that they are completely sym- 
metric with respect to all n! permutation of the indices 

V k h n = h n where V k € S n . (17) 

The operator Vk is the permutation operator that be- 
longs to the complete symmetric group S n . An important 
feature of this method is the availability of the analytical 
gradients of the total energy with respect to the geminal 
parameters. The gradients can be computed analytically 
and are given by the following expressions 



dg(i,2) 
db k 



dg(l,2) 
djk 



-b k r\ 2 e-^K 



(19) 



The AO integrals involving the gradients of the GTG 
functions are performed analytically and are computed 
with other AO integrals. 



B. Optimization of the trial wavefunction 

The optimization of the trial wavefunction is per- 
formed by performing a full configuration interaction 
(FCI) calculation on the congruent transformed Hamil- 
tonian. The FCI wavefunction is constructed by per- 
forming all possible excitations from the reference wave- 
function. [ ] This can be represented by the following 
expression, 

JVocc Wvir JVocc JV vil 

*fci = Co$+^^CX + ^^CK (20) 

a p a<b p<q 

N occ N vir 

+ E E 

a<b<cp<q<r 

where we have retained N v i T in the expression to em- 
phasize that only a finite number of terms are evalu- 
ated. This point will be a subject of discussion later 
in the derivation. The occupied and virtual orbitals 
are represented by (a, 6, c, . . . ) and (p, q, r, . . . ), respec- 
tively, and the CI coefficients are represented by (C£, . . . ) 
and are obtained variationally by minimizing the total 
energy. The construction of the full set of excitations 
and the determination of the CI coefficients are the two 
principle computational challenges associated with the 
FCI method. For very small molecules, the CI ma- 
trix can be explicitly constructed and diagonalized, how- 
ever, this simple approach becomes prohibitively expen- 
sive as the system size increases. Currently, there are 
various computational techniques for efficient calcula- 
tion of the expansion coefficients. [12, 46-50] The calcu- 
lation requires matrix elements involving the operators 
{(®k\O a \$ k >),a = 2, ... ,6} which are derived below. 

The matrix elements involving the 2-particle operators 
are evaluated as 

2! N occ 



(^o\0 2 \^o) = ^f^i-irHhiMPki^), (21) 

— 1 i x %2 
2! N occ 

<*o|02|*5) = ££(-ir<a*ilW^i), (22) 

A;— 1 ix 

(WW) = f>ir<a&|/i 2 |P feOT ). (23) 
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The matrix elements involving the 3-particle operators 



(18) 
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are evaluated as 



3! N, 



$o|0 3 |$o) = ^E E (-^ Pk {hi2h\h 3 \P k iii 2 k) 



k = l ■ 



3! JV D 



(24) 



*0|O 3 |*S) = ^£(-1)^^21/13^1*2), 



fc— 1 ii«2 



3! N Q , 



(25) 



(*o|0 3 |^)=^^(-ir(aW x |/ l3 |P fc Wi), (26) 

fc — 1 2l 

3! 

(*o|Os|*SD = £(-l) p *(a6cNi^r). (27) 



fc=i 



The matrix elements involving the 4-particle operators 
are evaluated as 



4! jv occ 

($oio 4 |$ ) = ^E E ( 28 ) 

i 21222324 

( — l) Pfe (?li2*3i4|^4|-Ffcil*2i3«4), 
, 4! Af occ 



fc— 1 HZ2Z3 

4! N occ 



(29) 



'2; 



k — l iii^ 



4! iV occ 



($o|04|*X) = E E(- 1 ) P M«^ 1 |^|P fc p?ni) ) 



fc — 1 2i 



(^o|0 4 |$E) - £)(-ir(aM^lflkP9»-*>. 



(30) 



(31) 
(32) 
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The matrix elements involving the 5-particle operators 
are evaluated as 



1 5! Nocc 

<4>o|0 5 |<I>o) = ^]r £ (33) 

/c=l 2122232425 

(-l) Pt (*'l*2*3*4*5|fr5|-Pfc*l*2*3*4*5), 
-, 5! W cc 

(*o|0 B |*S> = S E X! ( 34 ) 

fc=l 21222324 
(-l) Pfc (^12*3*4^5 |-PfcP*i*2*3*4), 

5! Af OC c 



(*o|0 6 |*S) = ^E E (- 1 ) Ph (o>biii2i 3 \h 5 \P kP qhi2i3) 

k—l 21*223 

(35) 

5! jV occ 
fc= 1 2^22 

(36) 



5! JV occ 

<*o|0 B |<TO = E E(-ir<atediil^lftP3r«i>, 
fe=i ii 

(37) 

5! 

(*o|o 5 |$s!) = E< aMe N p ^ s *>- ( 38 ) 
fc=i 

The matrix elements involving the 6-particle operators 
are evaluated as 

6! JVccc 

<*o|0 6 |<i>o) = ^£ £ ( 39 ) 

k=l 212223242526 

(-l) Pfc (*l*2*3*4*5*6l^6l^ , fe*l*2*3*4*5*6), 
, 6! Af OC c 

(*o|0 6 |$ p ) = ^£ Y, ( 4 °) 

fe=l 2122232425 

(-l) Pfc (mii2*3*4*5|^6|-PfeP*l*2*3*4*5), 
, 6! JV occ 

<*oiobi»3> = s E E ( 41 ) 

k—l 21222324 

(-l) Pfc (a6iii2*3i4|/i6|-PfeW*l*2*3*4), 
1 6! Af occ 

<*o|Oa|*£> = ^EE ( 42 ) 

fc=l 212223 

(-l) Pfc (a5ciii2* 3 |/i6|f , fcW* , *i*2*3), 

1 6! JV„cc 

<*o|06|*33> = ^£E («) 

ii2 2 

(-l) Pfc (a6cd*i* 2 |/i6|-PfeM* , 'S*i*2), 

6! JV occ 

(*o|0 6 |$ p ^> = E E (abcdehlhelPmrsth), (44) 

k—l %\ 
6! 

(*o|0 6 |*Kf> = Y,{^cdef\h 6 \P m rstu). (45) 

The computation of matrix elements in the above ex- 
pression requires atomic orbital integrals involving the 
GTG functions. One of the advantages of using the 
GTG functions is that all the AO integrals needed for 
the CT Hamiltonian calculation can be computed ana- 
lytically. Boys and Singer have derived the integrals in- 
volving GTG functions with s-type GTOs. Persson and 
Taylor have extended the method for higher angular mo- 
mentum by using the Hermite Gaussian expansion ap- 
proach. [43] Recently, Hofener and coworkers have also 
derived the geminal integrals by extending the Obara- 
Saika techniques for calculating the GTG integrals. [51] 
The solution for the CI coefficients requires diagonal- 
ization of the CI Hamiltonian matrix. However, the low- 
est eigenvalue and eigenfunction can be obtained without 
explicit construction and storage of the CI matrix. There 
are various efficient methods such as the Davidson diago- 
nalization to perform this task. [52] Recently, Alavi et al. 
have developed the FCIQMC method which allows very 
efficient evaluation of the FCI wavefunction. [53-57] 
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In the present calculation, the FCI eigenvector was ob- 
tained by performing the Nesbet update scheme and was 
selected because of its ease of implementation. [58] In the 
Nesbet method, a expansion coefficient is updated by 
Ac,,, 



where the update is calculated as 
Ac, = . °> , 



The energy is updated at each step using 

. „ cr„Ac„ 
AE = 



D + AD ' 



AD = Ac, 



2 ^ S^iCi + >S„ Ac, 



(46) 

(47) 
(48) 

(49) 
(50) 



The FCI energy can be recovered from the CT calculation 
by setting G = 1 



EpOi = lim i?cT- 

G — y 1 



(51) 



From the above relationship, we expect that the CTH 
energy calculated with G m i n will be lower than the FCI 
results. In the following section, we perform CTH calcu- 
lations on well-studied two-electron systems and compare 
calculated energies with reported benchmark values. 
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FIG. 1. Comparison of the exact ground state energy of 
Hooke's atom with the results from CTH and FCI calcula- 
tions. 



TABLE I. Geminal parameters for Hooke's atom using the 
6-3 11G basis set 



Number 


b k 


7fe 


1 


1.0000 


0.0000 


2 


-0.6090 


0.1050 


3 


-0.0709 


2.350 


4 


0.0216 


0.175 


5 


-0.0132 


1.120 



III. CALCULATIONS AND RESULTS OF 
BENCHMARK SYSTEMS 

The Hooke's atom is one of the few correlated two- 
electron systems for which the Schrodinger equation can 
be solved analytically. This feature has made it a test- 
ing ground for a wide variety of methods. [59-62] The 
Hooke's atom consists of two electrons in a parabolic po- 
tential. The Hamiltonian of that system can be written 
as 



H 



- \v\ + \kr\ + \kr\ + — 
2 2 2 2 7*12 



(52) 



where, all the quantities are expressed in the atomic 
units. The interaction between an electron and the nu- 
cleus is described with the harmonic potential. For k — 
0.25 a. u., the Schrodinger equation can be solved exactly 
and the ground state energy is equal to 2.0 Hartrees. [63] 
The Hooke's atom provides an ideal ground for testing 
the CTH method. The CTH calculations were performed 
using the 6-31 1G basis and the geminal parameters were 
obtained variationally from the solution of Eq. (4) . The 
energy was converged with respect to number of geminal 
parameters N s , and the results are presented in Fig. 1. 



It is seen that the energy was converged after addition 
of four geminal parameters and the optimized geminal 
parameters are listed in Table I. Comparing the energy 
with the exact result of 2.0 Hartrees, it is seen that the 
E[G min ] is slightly higher by 0.770 mHartrees (or 0.483 
kcal/mol). The optimized Slater determinant $ obtained 
in the previous step is used as the reference wavefunction 
for the CTH calculations and the results are summarized 
in Fig. 1. For G = 1, the CTH energy is identical to 
the FCI energy. However, inclusion of additional gemi- 
nal terms makes the CTH energy lower than the energy 
from the FCI calculation. It is seen that the CTH energy 
is in good agreement with the exact analytical results 
and is higher by 0.000296 Hartrees, these results are pro- 
vided in Table II. The CTH calculations were also carried 
out for the helium atom and the results are presented in 



TABLE II. Difference between exact and calculated energy 
Hooke's atom using the CTH method 

Hartree kcal/mol kj/mol eV cm" 
0.000296 0.186 0.777 0.00805 65.0 
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FIG. 2. Effect of basis set on ground state energy of Helium 
for HF, FCI, and CTH methods. 



TABLE III. Geminal parameters for Helium atom using 6- 
311G 



Number 


b k 


Ik 


1 


1.000000 


0.00000 


2 


-0.320260 0.57816 


3 


-0.063365 


10.3760 


4 


0.020918 


0.83536 


5 


-0.029282 


0.08799 



Fig. 2. The calculations were performed using different 
basis functions, and the results were compared with HF 
and FCI values. It is seen that for small basis sets, the 
E energy is lower than the FCI energy. We expect this 
because of the inclusion of the optimized geminal terms. 
The key result from Fig. 2 is that for small basis sets, 
the CTH method provides a substantial lowering of en- 
ergy with respect to the corresponding FCI values. The 
CTH calculations with respect to a small basis provides 
a wavcfunction that is comparable to the FCI wavefunc- 
tion at much larger basis functions. Since the cost of 
the FCI expansion increases sharply with the size of the 
underlying 1-particle basis, the CTH method provides an 
appealing alternative for obtaining accurate results when 
an FCI calculation is prohibitively expensive. The de- 
pendence of the CTH energies on the number of geminal 
parameters is shown in Fig. 3 and the optimized geminal 
parameters for the helium atom are listed in Table III. 

IV. DISCUSSION AND CONCLUSIONS 

The first geminal parameter is always set to b\ = 1 
and 71 = and is never optimized during the calcula- 
tions. When all the other N R — 1 geminal parameters 
are set to zero, these values of b\ and 71 represent the 



FIG. 3. Convergence of the CT Hamiltonian energy of the 
Helium atom with respect to the number of geminal functions. 
The calculations were performed with 6-3 11G basis set. 



G = 1 limit. Geminal parameters from &2 • • ■ &j\r K and 
72 • ■ • 7-/Vg are optimized to obtain G m - m as described in 
Eq. (4). This procedure ensures that the optimized en- 
ergy is always bounded from above by the HF energy. 
Figures 1 and 2 show the effect of inclusion of additional 
geminal parameters and it is seen that the second gem- 
inal parameter lowers the energy significantly. This is 
an important result and clearly indicates the importance 
of the geminal function in construction of the congruent 
transformed Hamiltonian. The set of {bk} was optimized 
without any constraint and it is seen from Tables I and 
III that the overal geminal parameter is negative. This 
is an expected result and is in agreement with previous 
work on explicitly correlated methods. [51, 64, 65] The 
negative values of geminal parameters indicate the role 
of the geminal function in providing a better description 
of the Coulomb hole. 

The analytical forms of the GTG functions are inher- 
ently approximate and are not capable of describing the 
cusp correctly because their first derivative vanishes in 
the limit of r ce = 

To assess the quality of the CTH energy, it is important 
to estimate how much of an error this feature introduces 
in the calculated energy. For the Hooke's atom this can 
be done in a a straightforward manner since the analyti- 
cal solution of the Schrodinger equation is known. From 
Table II, it is seen that the CTH energy is close to the ex- 
act ground state energy and is higher by 0.296 mHartrees 
or 0.186 kcal/mol. This difference between the CTH and 
the exact energy represents the upper bound in the er- 
ror that one can expect for this system by approximating 
the cusp with GTG functions. For the helium atom, the 
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TABLE IV. Comparison of ground state energy (in Hartrees) 
of the helium atom 



E 


function 


ref 


-2.900233 


FCI 


[75] 


-2.903041 


CTH 


this work 


-2.903724 


free ICI 


[69, 70] 



situation is less straightforward because we do not have 
access to the exact solution. Instead, we compared the 
CTH energies with other high-level methods from previ- 
ous studies [66-74] that include the exact cusp condition 
in the wavefunction. In order to achieve the best CTH 
energy, the calculation was performed with an aug-cc- 
pVTZ basis set and geminal parameters were optimized 
with respect to the aug-cc-pVTZ basis. Comparing the 
CTH method with the highly accurate ICI method by 
Nakatsuji [69, 70], it is seen that CTH energy is higher 
than the ICI energy by 0.429 kcal/mol. The compari- 
son of the CTH calculation to the ICI method and other 
highly accurate results can be seen in Table IV. The im- 
pact of electron-electron cusp on ground state energy was 
investigated in detail by Prendergast [12] and coworker 
using CI and QMC methods. Their study concluded that 
one can still expect to get mHartree level of accuracy 



even in situations where the exact cusp condition is not 
satisfied. Our study using GTG functions also confirms 
this observation. The use of GTG functions in the CTH 
method, represents a trade-off between the implementa- 
tion of exact cusp condition and analytical expression for 
computing the Gaussian-type geminal integrals. 

In conclusion, the congruent transformation of the elec- 
tronic Hamiltonian using Gaussian type geminal function 
is presented as a general method for calculating accurate 
ground state energy. The form of the congruent trans- 
formed Hamiltonian can be systematically improved by 
using the geminal function. It was found that a small 
number of geminal functions are needed to converge the 
energy. Futhermore, addition of just one geminal param- 
eter results in a substantial improvement in the accuracy 
of the wave function. For a given finite basis set the 
CTH energy was found to be lower that the FCI calcula- 
tion on untransformed Hamiltonian. The results indicate 
that the congruent transformed Hamiltonian provides a 
viable alternative for obtaining FCI quality energy using 
a smaller underlying 1-particle basis set. 
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